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^ ABSTRACT 

>■ We consider a three-dimensional bipolar force-free magnetic field with non 

^ zero magnetic helicity, occupying a half-space, and study the problem of its evo- 

^ lution driven by an imposed photospheric flux decrease. For this specific setting 

of the Flux Cancellation Model describing coronal mass ejections occuring in ac- 
(^ tive regions, we address the issues of the physical meaning of flux decrease, of the 

'^ influence on field evolution of the size of the domain over which this decrease is 
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imposed, and of the existence of an energetic criterion characterizing the possible 
onset of disruption of the configuration. We show that: (1) The imposed flux 
disappearance can be interpreted in terms of transport of positive and negative 
fluxes towards the inversion line, where they get annihilated. (2) For the par- 
ticular case actually computed, in which the initial state is quite sheared, the 
formation of a twisted flux rope and the subsequent global disruption of the con- 
figuration are obtained when the flux has decreased by only a modest amount over 
a limited part of the whole active region. (3) The disruption is produced when 
the magnetic energy becomes of the order of the decreasing energy of a semi-open 
field, and then before reaching the energy of the associated fully open field. This 
suggests that the mechanism leading to the disruption is nonequilibrium as in 
the case where flux is imposed to decrease over the whole region. 



Subject headings: magnetohydrodynamics (MHD) — stars: coronae — stars: 
magnetic field — stars: flare — Sun: coronal mass ejections (CMEs) — Sun: 
flares 



1. INTRODUCTION 

Magnetic flux cancellation (FC) is a well documented solar photospheric phenomenon 
(Welsch 2006) which plays the key role in the Flux Cancellation Model (FCM). FCM has 
been introduced initially to explain flux rope and prominence formation (van Ballegooijen 
& Martens 1989), and applied later on for modelling large scale eruptive events such as 
Coronal Mass Ejections (CMEs) occurring in active regions (see, e.g., Forbes et al. (2006) 
for a review). In the latter context, it has been actually developed in two settings differing 
from each other essentially by the way FC is introduced: FC is either taken to be a physical 
consequence of the photospheric turbulent diffusion, or just imposed as a time dependent 
boundary condition describing a decrease of the flux according to some prescribed law. In 
the former case, flrst considered in Amari et al. (2003b) and later on in Yeates & Mackay 
(2009), Yeates et al. (2010), and Aulanier et al. (2010), the dispersive effect of the turbulence 
on the flux of an active region leads to the bringing together on the inversion line of some 
amounts of flux of both polarities, which thus get annihilated by small scale mixing. This 
type of FCM is of course particularly relevant when one wants to explain the CMEs which 
occur during the decaying phase of an active region. The second type of FC implementation 
was flrst applied to study the possible disruption of 2D bipolar conflgurations (Forbes & 
Priest 1995), and later on of 3D bipolar conflgurations (Amari et al. 2000; Linker et al. 2001) 
and 3D quadrupolar ones (Amari et al. 2007). 

Although previous works have led to a good understanding of FCM in the "flux decreas- 
ing setting", there are many points which are still unclear. (1) The disappearance of flux 
on the boundary is imposed as a mere mathematical boundary condition, and doubts have 
been sometimes expressed about the possibility that it be related to some physical process. 
To keep viable this setting of the FCM, we thus need to provide a physical interpretation 
of flux decrease compatible with the observations. For instance, the latter show FC to be 
often associated to the mutual annihilation of opposite polarity flux elements brought into 
contact by photospheric motions, and an interpretation in terms of flux transport and anni- 
hilation on an inversion line would be certainly adequate. An other possibility would be to 
interpret FC in terms of emergence through the photosphere of either an U-loop or a bipolar 
loop, flux decrease occuring in the latter case once the magnetic axis has started emerging 
(see Fan (2001) and Amari et al. (2005) for simulations, and Lopez Fuentes et al. (2000) 



for observations). (2) Previously, flux decrease has been imposed to occur over the whole 
active region, or at least the whole central bipolar part in the case of a quadrupolar region 
(Amari et al. 2007). In many cases, however, one may expect FC to occur over only a small 
part of the region, and the question arises of the possibility of still triggering a disruption in 
that case. (3) If a disruption actually occurs, we have finally to address the question of the 
nature of its trigger. Basically, there are three possible mechanisms (Amari & Aly 2009): 
Non Equilibrium (NE), Quasi Non Equilibrium (QNE), meaning equilibrium too far to be 
accessible, and Unstable Equilibrium (UE). For instance, we found QNE to be responsible for 
the transition to very fast expansion exhibited by a flux rope twisted by boundary motions 
(Amari et al. 1996). In the case of a bipolar (Amari et al. 2000) or quadrupolar (Amari 
et al. 2007) configuration submitted to FC, on the contrary, we found NE to be at the origin 
of the eruption. This conclusion was established by showing that the magnetic energy of 
the system exceeds the energy of the associated totally open field in the former case, and 
the energy of a partially open field (having its open lines being connected only to the part 
of the boundary where strong currents once developed) in the latter one. Then we should 
look in particular for the existence of a similar energy criterion if a disruption occurs as a 
consequence of partial FC. 

The aim of the present Letter is to address these important issues in the case where the 
initial configuration is force-free and has a nonzero magnetic helicity. Physically, such a field 
may be thought of either resulting from the emergence of a twisted fiux rope (TFR) through 
the photosphere, or being the remnant of a configuration which has previously generated 
an eruptive event without fully relaxing to a potential field (Amari & Luciani 1999), or 
to be produced by photospheric twisting motions. It is worth insisting on the fact that 
the evolutions we consider are only driven by FC, in opposition to some other works (van 
Ballegooijen et al. 2000; Aulanier et al. 2010) in which FC is applied along with shearing 
motions. 



2. MODEL AND INITIAL CONFIGURATION 

In our numerical model, the "coronal half-space" {z > 0} above an active region is 
represented by a large computational box Qh = [—20, 20] x [—20, 20] x [0, 40], equipped with 
a non uniform mesh of 141x131x121 nodes. The used physical quantities are dimensionless. 
In particular the Alfven crossing time ta is taken as the time unit, flh contains a perfectly 
conducting low density plasma with an embedded magnetic field B. For t > 0, this system 
is brought into a two-stages MHD evolution: a preliminary stage, not claimed to represent 
an actual coronal evolution, which produces a force-free equilibrium with a nonzero helicity. 



and a second stage, which is the physically relevant one, in which the previous field is taken 
as the initial state of a FC driven evolution. Both phases are controlled by imposing on 
the "photospheric" boundary Sh of Qh (on which z = 0) the tangential component of the 
electric field, E^, a procedure which is known from our previous works to lead to a well posed 
problem. Quite generally, E^ can be Helmholtz decomposed according to 

cE, = V,0 + V,V^ X z, (1) 

where V^ = itd^ + 'y dy, and 0(a;, y, t) and ip{x, y, t) are solutions, respectively, of the equa- 
tions Vs^0 = cVs ■ Es and — Vs^-?/' = cz ■ V^ x E^. Using Faraday law thus leads to 

dtB, = -cV ■ (E, X z) = V,V- (2) 

The MHD equations are solved by our numerical algorithm (Amari et al. 1999, 1996). Small 
values are choosen for the dissipation coefficients: u = 10^^ — 10"'^ for the kinematic viscosity, 
and f] not larger than 10~^ for the resistivity, giving for our mesh resolution a Lundquist 
number not smaller than 10^. 

In the first stage, we start at t = from a bipolar potential magnetic field B^ = V\4-- 
Vj^ is the solution of a Dirichlet-Neuman boundary value problem for Laplace equation, in 
which we impose in particular the condition 

B^,{x, y, 0) = g(x, y) = e""'/"' L-^y-y^^/'^i - e-^y+y^^/'^'y] (3) 

on Sh, with yc = —0.8, (Tx = l,(Ty = 2. For < t < 400, we apply on Sh, as in Amari et al. 
(2003a), the tangential electric field Eg associated to ip = and a such that (p'iBz) = 
Bz^'[Bz), with ^{Bz) a prescribed function. This generates slow shearing/twisting Bz- 
preserving ideal motions of the magnetic footpoints at the velocity Vg = c'Vs<f) x z/i?-, 
(with m.&x{vs) = 10^^). A neighbouring nonlinear force-free equilibrium is next reached by 
performing for 400 < t < 800 = to a viscous relaxation phase during which (p = ip = 
and Bz is still conserved on Sh- As shown on Figure 1, this equilibrium is sheared along 
the neutral line, and because of twist, it exhibits away from that line the presence of strong 
electric currents correlated with the typical sigmoidal structure. It is worth noticing that, for 
this equilibrium, we have supp{FFF) C supp{Bz), where supp{FFF) and supp{Bz) denote 
the regions where, respectively, electric currents and Bz are strong enough. 



3. PARTIAL FLUX CANCELLATION 

For driving the FC phase, we fix E^ on Sh by setting (f){x, y,t) =0 and by taking ip 
to satisfy V/^(x,|/,t) = -fiC{y)Bz{x,y,0,to) = -fiC{y)q{x,y) on Sh, with /i = 10"^ > 0, 



and ({y) = [1 + tanh[(t/i — t/)/(i]][l + tanh[(t/ — i/q) /d\]/4:. ( exhibits a "plateau" of height 1 
in the interval [l/o?yi] and falls down to zero in a layer of thickness d. It is used as a mask 
controlling the size of the cancellation support, supp{FC). We use here the particular values 
yo = —0.5,yi = 0.5, d = 0.1, whence supp{FC) C supp{FFF) C supp{Bz) as the width 
of each Gaussian polarity is approximatly equal to 2 and decreases at a slower rate. The 
case considered in Amari et al. (2000) (where FC is enforced over the whole bipolar region) 
corresponds to C = 1- Using Eq. (2), we have dtBz{x,y,0,t) = — Ai^(l/)5'(a;,|/), whence after 
an immediate integration 

B,{x,y,0,t) = q{x,y)[l-fiay){t-to)]. (4) 

Then Bz suffers a linear decrease on Sh- But this does not lead to any specific problem as 
the reduction factor [1 — fi({y)(t — to)] > 0.62 > 0, the FC phase being limited to a duration 
of about 38 (see below). All along that phase, we regularly apply viscous relaxation runs 
to check if a neighbouring equilibrium exists or not. We found this strategy to be more 
appropriate than the one used in Amari et al. (2007), where the need to effect episodic 
relaxations was bypassed by choosing a much smaller value of /i (/i = 10"'^). 

As stated in Section 1, it is crucial for the validity of the model that the mathematically 
imposed flux decrease be interpretable in physical terms. Let us show that it can be actually 
considered to result from the transport of opposite polarities fluxes towards the inversion 
line / where they annihilate indeed. In any case, we can introduce on Sh an horizontal 
velocity of magnetic flux transport, u, by setting cE^ + B^u x z = 0. Using the Helmholtz 
decomposition (1) of cE^ and Eq. (4), we obtain at once 

u(a;, y, t) = -V,^(x, y)/[q{x, y){l - fiC{y){t - to))]- (5) 

During the FC phase, u appears to be just continuously rescaled by a time dependent factor 
on most of supp{FC), where ( = I, and to remain invariant outside supp{FC), where C = 0- 
Interestingly, we see that the cancellation flow support, supp{CF), coincides approximately 
with supp{Bz), and is then not reduced to supp{FC) for a localized FC. The velocity u 
computed at some time is shown on Figure 2a, and it is clearly seen to be directed indeed 
towards / in a neighbourhood of that line. Note that we have removed from the plot the 
parts very close to / and very far from the spots, respectively, where u becomes very large. 
The flux BzVl of B^ keeps however a flnite value on /, on which Vs"?/' does not vanish, and 
the flow thus continuously advects magnetic flux onto that line, where it disappears. On 
the other hand, it should be noted that the velocity pattern at any other time during the 
FC phase is identical to the one shown on Figure 2a because of the simple scaling property 
proved above. 

That flux decrease may be associated with flux annihilation on the inversion line / can 
be also simply proved in a quite general way as follows. Let us consider for a little while an 



arbitrary initial distribution Bzito) on Sh, and prescribe a flux decrease by setting dtB^ = 
fc{x, y, t), with /c > on S^, where -B^ > (/c < on S*^, where B^ < 0), and f^ feds = 0. 
We define as above the fiux transport velocity u by BzU = cEg x z = Vs0 x z — Vs'?/'. Then 
we have 

[B.UMdl = f{di(t) - dNij)dl = - f fjs > 0, (6) 

Ji Ji Js+ 

where N is the exterior normal on Sh to dS^ (and then to /), and un = u-N. For estimating 
the second member, we have used Gauss theorem, the equation V^^-?/' = fc (which results 
from Eq. (2)), and assumed d^ip and to vanish on the far part of the boundary. We can 
thus conclude that, "in the mean sense ", we have in a neighbourhood of / (into which N 
has been extended) un > on the + side and un < on the — side, which corresponds 
indeed to a finite amount of fiuxes of both signs being transported towards /, where they 
cancel. But of course, for arbitrary functions Bzito) and /c, and associated sinuous shape of 
/, we expect to have in general un < Q {un > 0) on the + (— ) side of some "httle" part of 
J, i.e., local emergence of fiux, and FC holds only globally. That this does not happen in the 
particular case considered in this Letter is due to the fact that we consider a configuration 
with antisymmetric fiux distribution {q{x, —y) = —q{x, y), whence I = {y = 0}), and impose 
fc = —f^Cq, with ( symmetric {C{~y) = C{y))- Finally, we note the following relation, which 
holds for the latter specific choice of fc- Consider an arbitrary curve C surrounding the 
positive part of supp{FC), with external normal N. Setting $(t) = J^^ ,p,^.^Bzds and 
proceeding as above, we get 

/ BzU^dl = -(ii$ = At / (qds ~ /i$(to) (7) 

Jc J supp(FC)+ 

Increasing the size of supp{FC) increases $(to) and thus the mean velocity < un >■ 

For a turbulent diffusion driven evolution (Amari et al. 2003b), the arguments above, 
based on our Helmholtz decomposition of Eg, also apply, thus providing a formal proof that 
cancellation fiows are associated to that process too. 



4. TWISTED FLUX ROPE, DISRUPTION, AND TRIGGERING 

MECHANISM 

We now describe the evolution of the coronal field driven by FC. 

(1) During a first phase, there is a transition from an arcade topology to a twisted 
fiux rope (TFR) topology. This transition is not instantaneous. Rather we observe more 



and more field lines close to / forming dips and a TFR progressively grows from those as 
flux is advected towards /. A coherent TFR exists at t = tpR- In particular J-shape loops 
have "merged", generating the inverse topology characteristic of the TFR, with dips able 
to support prominence material, as seen on Fig. 3. Such a feature was also found in the 
simulations reported in Amari et al. (2000). During the transition, the mutual helicity of 
the two J-shape loops is converted into the TFR self-helicity. 

(2) During this phase, viscous relaxation always leads to a neighbouring equilibrium. 
The TFR created along the neutral line grows up both laterally and vertically, as new flux 
participates to its structure. As expected for a system in near equilibrium (Aly 1984, 1991; 
Sturrock 1991), we find that WT,{t) <W(t) < W„{t), where W{t) is the magnetic energy of 
the configuration, and W.„{t) and W^{t) are, respectively, the energies of the potential field 
B-^(t) and the open field Bo-(t) having the same flux distribution on Sh as B(t). W{t) and 
WT,{t) are observed to decrease slightly, while the ratio W{t)/WT,(t) increases (see Fig. 2b). 
This behaviour is expected since WT^{t) depends only on the photospheric flux distribution, 
decreasing by FC, while W{t) depends also of the volumic distribution of coronal currents. 
As VFff(t), W„{t) depends only on the photospheric flux distribution and decreases, but the 
inequality Wait) > W(t) stays satisfied as FC has not yet been very effective. 

(3) After a critical time tgd ~ 838, the configuration experiences a global disruption, 
as shown on Figure 4, and FC is switched off. Viscous relaxation does no longer lead to 
an equilibrium close to B(t), and the configuration evolves dynamically towards opening 
(see Fig. 5). As in Amari et al. (1996, 2000, 2003a,b, 2007), opening is characterized by 
a transition to very fast expansion suffered by a bundle of lines, which thus close down 
eventually at very large distances. Along with that process, reconnection develops through 
the overlying arcade and a current sheet forms below, associated to dissipation. At this 
stage, since our results rely only on viscous relaxation, it can be either that there exists an 
accessible equilibrium, but that it is unstable (case UE of our classification) or that such an 
equilibrium does not exist (case NE of our classification). Discriminating between the two 
situations could be done in principle by first trying to compute force-free equilibria satisfying 
the boundary conditions on i?„ and a = (V x B)z/Bz given by the evolution model before 
and after tgd, respectively (e.g., by using the Grad- Rubin methods developed in Amari et al. 
(2006)), and by studying their stability when they exist. But this is a difficult task since 
convergence of the equilibrium algorithms might be difficult to achieve around the critical 
point, while testing stability requires superimposing perturbations of some types. Then, 
as we favor a priori NE to be responsible for the disruption, we adopt as in Amari et al. 
(2000, 2007) the strategy consisting to show that a necessary condition for the existence of 
an equilibrium ceases to be satisfied at tgd- 



(4) The idea is that it becomes energetically favorable for the system to open once the 
magnetic energy W{t) starts exceeding the energy of an accessible semi-open field Bso{t)- 
The Bso{t) which seems to be relevant here is defined as follows: it has all its lines connected 
to supp{SO) being open and all the other lines being closed. A priori supp{SO) can be 
reasonably defined in two different ways: it can be identified either with the part of Sh to 
which the fast expanding lines are connected, or with the part of Sh where the twist/shear 
of the initial force-free field is concentrated. Happily, however, the two definitions select the 
same region supp{SO). We thus estimate Wsoif) by equating it to the energy of the potential 
field B'^(t) satisfying B'^^ = |go| on supp{SO) and B'^^ = go on Sh\supp{SO). As guessed, it 
appears that the transition to NE occurs indeed when Wso{t), which decreases faster than 
W{t), becomes comparable to the latter. In fact Figure 2b shows that, while W^(t)/W^^(t) 
increases, Wso(t)/WT,(t) decreases, and the two curves cross for t ~ tg^. From this point, 
any attempt to find a neighbouring equilibrium fails and the configuration experiences a 
major disruption as shown on Figure 5. 

We note finally that the amount of cancelled flux at tg^ is about 6%, a small value 
compared to the 30% found in Amari et al. (2000) and even below the 9% found in the 
quadrupolar case in Amari et al. (2007). Note that the initial states considered in these 
previous works were also quite highly sheared. We can thus state that a small amount of FC 
is sufficient to trigger a large disruption. But we do not pretend in any case that the 6% value 
is general. Depending on the initial state, a higher or lower amount of partial cancellation 
may be needed for a large disruption to be produced. And for low shear initial states, it is 
even possible that only a confined eruption be obtained as in Amari et al. (2003b). 
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Earth Connection Theory Program, NASA's STEREO/SECCHI Consortium and Centre 
National d'Etudes Spatiales. The numerical simulations performed in this paper have been 
done on the NEC SX8 supercomputer of the Institute I.D.R.I.S of the Centre National de la 
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Fig. 1. — Selected field lines of the initial force-free configuration reached after a shearing- 
twisting phase followed by a viscous relaxation. Strong shear is accumulated along the 
neutral line. 
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Fig. 2. — (a) Plot of the velocity field associated with the imposed flux decrease. This flow 
is indeed a cancellation flow as it brings flux elements of opposite polarities into contact 
along the inversion line, where they get annihilated, (b) Time variations of some important 
quantities during the phase of twisting by boundary motions of the initial potential field 
(0 < t < 400), the phase of viscous relaxation (400 < t < 800), the phase of FC (800 < 
t < 838), and the nonequilibrium phase (838 < t). Are represented the normalised unsigned 
flux F{t)/FQ through Sh, the free magnetic energy measure Vr(t)/PVV(t) (which increases 
during the FC phase), and the free energy measure Wsoit)/W.^{t) (which decreases during 
FC phase) of a semi-open field which has all its open lines originating from the region on 
which the shear /twist is distributed. The last two curves intersect at some critical time, 
corresponding to a FC of less than 6%, beyond which the global disruption occurs. 
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Fig. 3. — Selected field lines of the equilibrium configuration obtained at t = 820, i.e., during 
the phase of localised flux cancellation. A TFR exhibiting dips favourable to the support of 
cool material has formed from the two J-shape loops seen on Figure 1. 
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Fig. 4. — Selected field lines of the evolving configuration at t = 838, further during the 
phase of localised flux cancellation . Beyond this state viscous relaxation to an equilibrium 
does no longer hold when flux cancellation is switched off. 
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(a) 



(b) 



Fig. 5. — Selection of field lines of the configurations obtained at (a) t = 885 and (b) 
t = 1135 from the configuration of Figure 4, with FC switched off. The global disruption 
involves opening, reconnection through the overlying arcade and below, and formation of a 
current sheet, associated to a high dissipation of magnetic energy. 



